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Abstract 



A construction of the bi-Hamiltonian structures for integrable systems on regular 
time scales is presented. The trace functional on an algebra of 5-pseudo-differential 

■ operators, valid on an arbitrary regular time scale, is introduced. The linear Poisson 
I tensors and the related Hamiltonians are derived. The quadratic Poisson tensors is given 

■ by the use of the recursion operators of the Lax hierarchies. The theory is illustrated by 
. ' A-differential counterparts of Ablowitz-Kaup-Newell-Segur and Kaup-Broer hierarchies. 

00 . 

O 1 Introduction 

> 



The concept of integrable systems on regular time scales can build bridges between field 



' systems and lattice systems. This concept provides us not only a unified approach to study 
^ - on discrete intervals with uniform step size (i.e. lattice Hli) and continuous intervals but also 
an extended approach to study on discrete intervals with non-uniform step size (for instance 
g-discrete numbers Kg) or combination of continuous and discrete intervals. 

The approach of time scales allows the unification of such classes of nonlinear evolution 
equations like field soliton systems [H [21 El IH E] , lattice soliton systems O El El E] , g-discrete 
soliton systems [lOl [HI [121 [13] and others. The above approach was initiated in [Hj where the 
Gelfand-Dickey construction was extended. The theory was further developed in [15], where 
systematic construction of (1 -|- l)-dimensional integrable systems on regular time scales was 
presented, and a very effective tool, classical i?-matrix formalism, was utilized. The i?-matrix 
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formalism provides a construction of infinite hierarchies of mutually commuting vector fields. 
In [15], we examined the general classes of admissible Lax operators, presented examples of 
integrable systems on time scales that can be written in an explicit form. We also explained 
the source of constraints first observed in [Tl] . 

The greatest advantage of the classical i?-matrix formalism is that it allows the construc- 
tion of the bi-Hamiltonian structures and conserved quantities. The goal of this work is to 
present bi-Hamiltonian structures for A-differential integrable systems on regular time scales. 
Thus the main result of this article is the formulation of an appropriate trace form on the 
algebra of (5-pseudo-differential operators, that is valid on an arbitrary regular time scale and 
in particular the real case recovers the trace form of pseudo-differential operators [2]. If the 
appropriate constraints are taken into consideration then the trace form recovers also the one 
of shift operators [7] . 

In Section [2l we give a brief review of the concept of time scales, including A-derivative 
and A-integrals. In Section [3l we fix the class of the A-differential evolution equations under 
consideration. Besides we define appropriate functionals and their variational derivatives. In 
SectionsHlandEl we describe the algebra of (5-pseudo-differential operators, the Lax hierarchies 
and the constraints that appear naturally between the dynamical fields of admissible finite- 
field Lax operators. In Section [6], in order to find the bi-Hamiltonian structures, we introduce 
a trace functional on the algebra of 5-pseudo-differential operators in terms of which we 
construct the linear Poisson tensors and the related Hamiltonians. The quadratic Poisson 
tensors are reconstructed in the frame of the recursion operators [T3] of the Lax hierarchies. 
Finally, in Section [TJ the theory is illustrated by bi-Hamiltonian formulation of finite-field 
integrable hierarchies on regular time scales which are A-differential counterparts of Ablowitz- 
Kaup-Newell-Segur (AKNS) and Kaup-Broer hierarchies. 

2 Calculus on time scales 

A time scale T is an arbitrary nonempty closed subset of real numbers M [T6| [T71 [T8| [T9] . For 
the definition of the derivative on time scales, we use forward and backward jump operators 
a, p : T ^ T defined by 

cr(x) = inf {y G T : y > x} p{x) = sup {y : y < x}. 

We set in addition cr(maxT) = maxT if there exists a finite maxT, and p(minT) = minT 
if there exists a finite minT. The jump operators a and p allow the classification of points 
on a time scale in the following way: x is called right dense, right scattered, left dense, left 
scattered, dense and isolated if a{x) = x, a{x) > x, p{x) = x, p{x) < x, a{x) = p{x) = x 
and p{x) < X < cr{x), respectively. Moreover, we define the graininess function /i : T ^ T as 
follows 



Besides, T*^ denotes a set consisting of T except for a possible left-scattered maximal point. 
Set = minT if there exists a finite minT, and set = —oo otherwise. Also set x* = maxT 
if there exists a finite maxT, and set x* = oo otherwise. 

Let / : T — ^ M be a function on a time scale T. Delta derivative of / at x G T'^, denoted 
by A/(x), is defined as 



cr(x) — X. 



A/(x) = lim 




s G T, 



provided that the hmit exists. A function on a time scale is said to be A-smooth if it is 
infinitely A-differentiable at all points from T**. 

If functions /, (7 : T ^ R are A-differentiable, then their product is also A-differentiable 
and the following Lebniz-like rule holds 

A(/^7)(x) = gix)Afix) + fiaix))Agix) 

= f{x)A9{x)+g{a{x))Af{x) ' 

Besides, if / is A-differentiable function, then 

/Kx)) = /(x)+/i(^)A/(x). (2) 

If a; G T is right-dense, then = and the relation ([2]) is trivial. 
The shift operator E is defined by the formula 

Ef{x) = xeT. 

Moreover the relation ([2j) implies that 

E = l + fiA. (3) 

In particular, if point x G T lies within some continuous interval, being part of a time 
scale, or if the time scale T = M, then A-derivative is ordinary derivative with respect to x, 
i.e. A = dx- If X e T is such that fi{x) 7^ 0, then A = j^{E — 1). This is the case when x is 
an isolated point, for instance T = Z or Kg. 

Every continuous function / : T ^ M possesses A-antiderivative F : T ^ M such that 
AF{x) = f{x) holds for all x G T'^. Thus we define A-integral from a to 6 of / by 

f f{x) Ax = F{b) - F{a) a,beT. (4) 

J a 

Notice that, for every continuous function / we have 

r''^^ fix) Ax = ii{x)f{x). (5) 

J X 

Hence, it is clear that the A-integral is determined by local properties of a time scale. 

In particular, when the points a and h lie within continuous interval, being part of a time 
scale, then (jlj) is an ordinary Riemann integral. If all the points between a and h are isolated, 
then h = a"' (a) for some n G and A-integral is a sum (this follows immediately from ([5])), 
i.e. 

I n— 1 

fix) Ax = 5^/xK(«))/(^^(a)). 

i=l 

For more complicated time scales which are combinations of continuous intervals, isolated 
points, etc., the integrals can be constructed by appropriate gluing of Riemann integrals and 
sums. 

The integration by parts formula follows from the Leibniz-like rule ([I]) as 

f f{x)Ag{x) Ax = f{x)g{x)t - f g{<j{x))Af{x) Ax, (6) 

J a J a 



where / and g are continuous functions. The generahzation of (jl]) to the improper integral is 
clear. Thus, we define A-integral over an whole time scale T by 



/ /(x) Ax := / /(x) Ax = lim F{x) — lim F{x) 

Jj J^^ x^x* x^x, 

provided that this integral converges, i.e. the limits exist. 

For our purposes, we demand the time scales where the forward jump operator a : T ^ T 
is invertible. A time scale T is called regular if a{p{x)) = x and p{a{x)) = x for all x G T. 
The first condition implies that a is 'onto' and the second condition implies that a is 'one- 
to-one'. Thus on a regular time scale a~^{x) = p(x). Actually, a time scale is regular if 
and only if each point of T \ {x*, x*} is either two-sided dense or two-sided scattered and the 
point X* = minT is right dense and the point x* = maxT is left-dense [l^. For instance 
T = [-1,0] U : A; G N} U {k/{k + 1) : A; G N} U [1,2] is a regular time scale. 

Let us consider some particular examples of regular time scales: 

The real case, T = M. We have cr(x) = x and n{x) = for all x G M. In this case 
A-derivative and A-integral are such that 



A/(x) = dj(x) and f /(x) Ax = / 



/(x) dx. 



The lattice case, T = H7j. Let hhe a positive parameter. In this case cr{x) = x + H and 
/i(x) = h, where x G ^Z. A-derivative and A-integral have the form 



Af{x) = Uf{x + n)-f{x)) and / /(x) Ax = ^V/(n^). 



nez 



The g-discrete numbers, T = := U {0} (g > 1). For x G Kg, one finds that 
cr(x) = gx and /i(x) = (g — l)x. Then 

A/(x)= ^(^^)-^(^) 



(g — l)x 
where x 7^ 0, and 

f /(x) Ax = 5^g"(g-l)/(g"). 

3 A-differential systems 

Consider A^-tuple u := (tii, . . . , Mat)""" of dynamical fields : T — > M being A-smooth functions 
on a regular time scale T. Let 

C = {Auk ■.k = l,...,N;AeS}, 

where 

S = |A*^At^'^ ■ . . . ■ A*"At^" : n G No,2i, ji . . .,tn,Jn G n} , 



and A"^ is defined by ( llOp . Therefore S is tlie set of all possible strings of A and A"^ operators. 
Note that A and A^ do not commute. 

What we mean by a A-differential system, is a system of evolution equations 

ut = K[u], (7) 

where t G M is an evolution parameter (time), Ut := ^ and K := {Ki, K2, . . .)""" with Ki being 
finite order polynomials of elements from C, with coefficients that might be time independent 
(A-smooth) functions. 

Additionally, we assume that all fields u with their A-derivatives are rapidly decaying 
functions as x goes to or x*. Then, the functionals have the following form 

F{u) = [ f[u] Ax, (8) 

where f[u] are polynomial functions of C. Clearly two densities give equivalent functionals 
([H]) if they differ modulo exact A-derivatives. Having defined the class of evolution systems 
(I7j) and functionals (E]), we further proceed in a standard way, that is we define the duality 
map, Poisson tensors, etc. 

The integration by parts formula (E]) leads us to the relation 

[ Aif)g Ax = - [ fAE-\g) Ax =: f f A\g) Ax. (9) 
Jt Jt Jt 

Thus the adjoint of A-derivative is given by 

A^ = -AE~\ (10) 

Note that 

= 1 + fiAl 

Besides, by the use of ([3]), one finds that 

{E^y = n{l + fxA)^ = /i - /iAE^V = f^-{E- 1)^~V = ^"V (11) 
Consequently, the variational derivative of a functional in the form ([H]) is defined by 

|^ = Ea'K^ *=1.....A'. (12) 

Suk ^ d{Auk) 

Notice that ^A = 0, therefore the definition of variational derivative (fT2l) is consistent with 
the definition of functionals (Ml ■ 



4 (5-pseudo-difFerential operators 

We introduce 6 operator acting on A-smooth functions u : T 

6u := An -|- Eu6, 



(13) 



which is consistent with the Leibniz-hke rule ([T]). By the use of fll3p . we have 

= E-^u6~^ + E'^A^u6~^ + E-^A^\6-^ + ... . 
The generahzed Leibniz rule for the 5-pseudo-differential operators takes the form 

oo 

sy = J2s'kf^'"' ^ez, (14) 

k=0 

where 

gn ^ ^k^n-k + _ , + ^n-fc^fc 72 ^ ^ 0, 

is a sum of all possible strings of length n, containing exactly k times A and n — k times E; 

5^ = E-^ (^A^''E''+^ + ... + E"+^A^''^ for ri < and A; ^ 

consists of the factor E~^ times the sum of all possible strings of length k — n — 1, containing 
exactly k times and —n — 1 times E~^; in all remaining cases S*^ = 0. Besides, we have 
the recurrence relations 

5^+1 = S'^E + ^^lA for n^O (15) 

and 

k 

Lemma 4.1 For all ?i G Z, the relation 

j2i-firs-, = {E-f^Ar = i (17) 

holds. 

The proof is postponed to Appendix. 

When X G T is a dense point, i.e. fi{x) = 0, then the rule ( fT4l) is of the form 



fc=0 ^ ^ 



^kj^n-k n^%, 



where (^) is a binomial coefficient such that (^) = "("~^)'--'^("~^+^) ^ ^^^^ particularly when x 
is inside of some interval then A = d^- Thus, in this case we recover the generalized Leibniz 
formula for pseudo-differential operators. 

For X G T such that /i(x) 7^ it is more convenient to deal with the operator ^ := /i5 
instead of 5. By the use of ( TT3i) . the generating rule yields 

iu= {E -l)u + Eu^, 

and hence it follows that 



^"^ = E it) (^ - i)'i?"-vr~' nez. 

k=0 ^ ^ 



(19) 



The important fact is that the operator A = aj(5* has a unique ^-representation A = 
a^^\ and nonnegative (negative) order terms with respect to 6 transform into nonnegative 
(negative) order terms with respect to C,- 



I 



5 Lax hierarchies 

We define the associative algebra of (5-pseudo-differential operators 



= d^k © 0<fc 



i^k ) K. i<k 



equipped witli a Lie bracket given by tlie commutator [A, B] = AB — BA, wliere A,B E Q. 
The classical /^-matrices following from the decomposition of g into Lie subalgebras are 

= 2*^-^^^ ~ ^<'=) ^ ^^i" ~ 2^ 2~ ^'^'^^ 
where = or 1 and projections are such that 

A^k = Y,<^^^' for A = 5^a,5\ (21) 

i^k i 

According to the classical i?-matrices fl20l) we have two Lax hierarchies of commuting 
evolution equations [15] 

{L^)^,,l\ A; = 0,1 nGN, (22) 

generated, in general, by fractional powers of some Lax operator L. The general admissible 
finite-field Lax operators have the form 



L = un6^ + un-i6^~^ + ... + U16 + U0 + + ^,5" Vs, (23) 

s 

where for k = the field mat is a nonzero time-independent field and U-i = 0. Analysing (l22l) 
one finds that 



for details see [15]. Hence, there arises natural constraint between dynamical fields from (1231) 
given by 

N+k-l 

J2 {-fif^'-'-'u. + i-nf-' Yl ^^^^ = (24) 

i=—k s 

where a is time-independent function (for k = 1 nonzero when /i = 0)0. Notice that, the 
constraint (IMl) is compatible with the dynamics of Lax hierarchies ( l22l) . Using (12^ one can 
eliminate one dynamical-field. The convenient choice is to eliminate the field u^^i for A; = 
and the field mtv for A; = 1. 

It is clear that in the case of T = M the Lax hierarchies (I22p yield field soliton systems 
and in particular for Lax operators in the form ( l23i) one recovers the results of [3], H] . In the 
case of T = Z one obtains lattice soliton systems that are equivalent to the ones considered 
in [9]. Notice that in [9] the i?- matrix for = 1 has a slight different form, this follows from 
the fact that the construction in [U] is by means of shift operators. Similarly for T = Kg the 
Lax hierarchies leads to g-discrete soliton systems, etc. 



^Notice that in the formulae of Theorem 3.5 from [15j there arc misprints in degrees of powers of fi. 



6 Hamiltonian structures 

Let A = aj5* be a (5-pseudo-differential operator. We define the trace form by 

TrA := - / - A^o\s=-i = [ J^i-fiy-'ai Ax, (25) 



i<0 



wliere A<o is tlie projection onto negative terms. To show that the substitution 5 = — ^ in 



( 12511 is well-posed, we state the following proposition. 

Proposition 6.1 Let the 5 -differential operators A and B be such that (y4i?)<o = AB , then 

1 f I 

- AB\^^_i Ax= - A\s^_i B\^^_i Ax. (26) 



Proof. It is enough to consider A = aS"^ and B = b6"' such that m + n < 0, thus 



Ax 



A/^ '""^ A/i 

= / a V(-/i)'=-™-"-^^^6 Ax = / a6(-/i)-™""-^ Ax, 

where the equality (fTTl) is used. Consequently ( !26l) follows. □ 

Roughly speaking, the above proposition implies that the multiplication operation in the 
algebra g of (5-pseudo-differential operators commutes with the substitution 5 = —- given in 
the trace form 



Note that, in particular for the points x G T such that /i(x) = 0, the trace form turns 
out to be 

TiA = / a_i Ax, 



T 



where A = aj5''. Thus when T = R, we recover the trace formula for the algebra of pseudo- 
differential operators [2]. For the case /i(x) 7^ 0, the trace form ( l25l) within the algebra of 
^-operators is given by 

TtA := - [- A<o|^_i Ax = - / - ^(-1)^< Ax, 

where A = J2i (^'iC- 

Theorem 6.2 The inner product on g defined by the bilinear map 

{■,■),: gxQ^K {A,B),:=Tt{AB), (27) 
in terms of the trace ( l25l) . is nondegenerate, symmetric and a.d-invariant, i.e. 

iA,[B,C])^+i[B,A],C)^ = 0. 



Proof. The nondegeneracy of (1271) follows immediately from the definition of the trace. 

In order to show that ( l25l) is symmetric, it is enough to consider the monomials A = a5™ 
and B = M". Then, if m, n ^ 0, we have Tt{AB) = Tt{BA) = 0. If m, n < the symmetricity 
immediately follows from (!26|) . Thus, it remains to prove the case when one of the operators 
A and B is of positive order and the other one is of negative order. 

Without loss of generality, let m > and n < 0. We consider the cases fi{x) = and 
7^ 0, separately. 

For n{x) = 0, we have 



Tt{AB) = Tr (a5™M") = Tr I V (™) aA'=M"*+"-^' 



.fc=0 



[ ( ^^JaA^+^+ife Ax. 

/ Vm+n+l/ 



The converse formula for (fTSl) has the form 



«<5- = ^5"-^Y^jAtV 

k=0 ^ ^ 



(2^ 



Hence 



X 



Ti{BA) = Tr(M"a5'") = Tr (™)M"'+"-^A^^a j = J {^^^^^)bA^''"'^"^'a A 
= / LZ+>^"'^''^'bAx = TT{AB), 



where we make use of (Q. 

For n{x) 7^ 0, we pass to the calculations in terms of ^-pseudo-differential operators. Let 
A = a^™' and B = b^^ with m > and n < 0. We have 



Tt{AB) = Tr (aC^D = Tr ( ^ ~ ifE^^-^i 



m+n—k 



a-0 



-/^ E ©(-1) 

'^'^ ^ k=m+n+l 



m+n—k 



a(E - ifE'^-H Ax 



The converse formula for ([T9l) is given by 



k=0 ^ ^ 



Let f{E) be a polynomial function of E. Then by (fTTl) . it follows that 



Therefore 

/ m 

1 ^ wk 



Tl{BA) = Tr ipCaC) =^^\Y^ ^rn-^J^^m+n-k ^^-1 _ 1)'^ ^fe-m^ 

\fc=0 

= - / - V (E-i - l)^'E'=-™a Ax 

= - [ - Yl (fc) - l)''^"'"''^ Ax = Ti{AB). 

The symmetricity of the trace functional on the algebra of ^-pseudo-differential opera- 
tors implies the symmetricity of the trace functional on the algebra of (5-pseudo-differential 
operators for fi{x) ^ 0. 

Hence, the inner product ( !27|) is symmetric. Finally, the adjoint invariance of (1271) follows 
from the fact that the inner product (1271) is symmetric and the multiplication operation de- 
fined on the algebra % of 5-pseudo-differential operators is associative. □ 

For the next proposition, we have to consider (1 + and its expansion such that it is 

valid for all points of T including the case /i = 0. Besides, we assume that the expansion of 
(1 + fi6)~^ is given by nonnegative order terms in 5-pseudo-differential operators. We derive 
the following expansion 



oo 



;i + := $^(-5)^/ + A/+^) = (29) 

/i 



k=0 k=0 



which can be verified by multiplying both sides of the expression ( 129|) with (1 + fi6) from 
right-hand side. Hence 



1 + ^5)-^(i + ^5) = + 

fc=0 ^ k=0 
k=0 ^ k=0 k=0 

= ^ - A/. + ±{-6^ - / - aM = 1, 



k=l 

where (El) and the converse formula (l28l) are used. 



Proposition 6.3 T/ie alternative formula of the trace form ( 1251) is given by 



TrA = f ^-^les {A{1 + fiS)-^) Ax, (30) 



where 

resA := a_i /or A = aj(5* 



Proof. We first calculate the residue. Thus 



res {A{1 + fi5)-^) = les IJ^ Zl^"^)''"^'^^ 



+k ' 



k=0 i ^ 



I + . . . j = res [J2{-ir~'a,^6~' + 



. i<0 / \i<0 



Substituting the residue into the alternative trace form ( 1301) . we obtain 

"es [-^[1. -\- fio) ) lAx = I 

Hence the trace forms fl^ and flHUjl are equivalent. □ 



TtA = [ ^^res {A{1 + /i^)"^) Ax = [ V(-/i)-*-^ai Ax. 



Notice that, the definition of the trace form given by ( l30i) is very similar to the trace 
formula introduced in [H]. However the trace form in [14] is valid only either on T = M or 
on regular-discrete time scales, i.e. on such regular time scales where all points are isolated 
and /X 7^ 0. One of the main contribution of this article is to generalize the trace form given 
in [13]. The trace form ( |25l) is valid for an arbitrary regular time scale. In particular when 
T = M, the trace form (l25l) implies the standard trace formula for the algebra of pseudo- 
differential operators. Thus, for the use of the alternative trace form (l30l) . we have to choose 
that (1 + fiS)~^ expands into nonnegative order 5-operators fl^^ and the expansion is also 
valid when = . 

Observe that, one can define alternatively the following trace form 

Ti'A := [ - A^o\s=-i ^x = - [ V(-/i)-^-ia, Ax, (31) 

valid on regular-discrete time scales, that is for 7^ 0. Choosing the expansion of (1 + fiS)~^ 
into negative order terms, 

00 _ 

(1 

and 7^ 0, the formula ( 1301) yields (l3Tll . i.e. one recovers from (l30l) the trace formula of shift 
operators, as in [H]. 

The shift operator can be introduced by the relation: S = 1 + fi5, then S"^u = E"^uS"^ for 
m e Z. The expansion of the operator A by means of shift operators S, i.e. A = a[S^ (we 
assume that expands into negative order terms of shift operator S) allows us to obtain 
from ( 13T|) the standard trace form of the algebra of shift operators 

Tr'^ := / -a' Ax. 

Jt 

Although the traces fl25l) and fl3T|) are not equivalent in general, they are closely related 
to each other on regular-discrete time scales. To be more precise, on regular-discrete time 



scales when applied to the constrained operators such that = const, the traces 0251) 

/J 

and (IHTj) are equal up to a constant, as 

^^ol5=„i = - A<o\s=^i + const. 

Notice that, this is the case of the Lax operators ( l23l) for which the constraints ( l24l) are taken 
into consideration. Thus it is clear that in the case of lattice time scale, i.e. T = Z, the 
Hamiltonians defined below, as well as the following scheme, are equivalent to the one from 
[0]. By similar observations one also finds that for T = Kg one recovers from fl30p and (1311) 
the trace form of g-discrete numbers, see the appendix of [12] . 

To sum up, the trace form (l25l) is valid on arbitrary regular time scales and in particular 
for T = R give the standard form of pseudo-differential operators. Besides, if the appropriate 
constraints are taken into consideration, fl^3|) also recovers the trace forms for T = Z of 
'lattice' shift operators and for T = Kg of g-discrete numbers. 

In order to define the Hamiltonian structures for ( l22l) we first need to find adjoint of R- 
matrices ( l20l) . i.e i?^, such that {A, RB)^ = {R^A, B)q. Using the alternative definition of the 
trace form ( 130|) one finds that 

R^ = pU-\ A; = 0,1, (32) 

where 

Pl^A={A{l+^^5r')^^^{l + ^^5) (33) 
and the projections are such that 

B<^k=Y.^^^^ for B = Y,^%. 

i<~k i 

Notice that the above projections are defined on the operators given in a different form than 
in dnD. 

The existence of the well-defined inner product ( 1271) allows us to identify g with its dual q*. 
Let JF(0 = g*) be the space of smooth function on q consisting of functional (E]). The linear 
Poisson tensor has the form [20l [5] 

TTodH = [RdH, L] + R^ [dH, L] 

= [L, dH^k] + {[dH, L] (1 + fxS)-') (1 + fxS) k = 0, 1, ^^^^ 

where H G J^{q). 

We do not present explicit form of the differentials dH with respect to Lax operators (l23l) 
as it would be cumbersome, but we explain how to construct them. We postulate that 

n 
i=l 

where n is the number of independent dynamical fields in ( 123|) (the number of the rest of the 
dynamical fields after taking the constraint ( l24l) into consideration). Thus, we look forward 
to express 7i's in terms of dynamical fields of (123|) and their variational derivatives by the use 
of the assumption 




Some useful formulae to calculate linear Poisson tensors ttq for the examples are derived 

as 

P^Q^aS^^b) = aS^^b + fiab 
Pl^{a6-^b) = a6-^b-6-^ab 
Pl^{6~^a6~^b) = 6~^a6~^b + 6~^fiab. 

The case of the quadratic Poisson tensor tti is more delicate and for the construction 
there appears additional conditions on R and -R^ [HI El [201 E]- The form of (IH^ with 
does not allow us to proceed in a standard way anymore, thus we omit the construction of 
the quadratic Poisson tensor from the i?-matrix scheme. In article [15], we have constructed 
recursion operators $ for the Lax hierarchies fl22l) . such that 

$Lf = Lt 

Thus, as we know ttq, the quadratic Poisson tensor vri can be reconstructed alternatively by 

The recursion operator $ is hereditary at least on the vector space spanned by the symmetries 
from the related Lax hierarchy. Therefore the Poisson tensors ttq and tti are compatible [2TII5]. 
Hence, the Lax hierarchies fl22l) have bi-Hamiltonian structure 



where the related Hamiltonians are given by 
They are such that dH„ = . 

7 Examples 

A-differential AKNS, A; = 0. The Lax operator (1231) for = 1, with the constraint 
(a = 0) is of the form 

L = 5 + fiipif + #^ V- (36) 
The first and the second flows from the Lax hierarchy (1221) are 

ipti = /i^'V + ^i'^ 

and 

For Lax operator ( !36l) the differential of an functional H, such that (j35l) is valid, is given 

(p dip ip \^ / tfxp + jjtpA'ip 



(37) 



where 

A-iP—- — 

and is a formal inverse of A. Then, one finds the hnear Poisson tensor (E 
The recursion operator constructed in [T3] has the form 



1 
-1 



A + 2/i?/'v9 + 2^A~V /i^^ + 2?/'A~V 
-liLp^ - 2ipA~^ip At - 2^A^ V 

Hence, the quadratic Poisson tensor is 

= = ^ -At + 2^A- V -/i^2-2^A-V J- ^ ^ 

The skew-symmetricity of (1391) follows from ([3]) and ( ITOl) . The first three Hamiltonians are 



Hq = I ipip Ax 
1 



T 

/ 1 

2,„2 



H-2 = j Q/i^V^V^ + V^V^ + '/'A^?/^ + fiipip'^Aip + /iV'VAV^ Ax 



Particularly, when T = M the above bi-Hamiltonian hierarchy is exactly the bi-Hamiltonian 
field soliton AKNS hierarchy [1]. In this case the first nontrivial flow is the second one ( |38l) . 
i.e. the AKNS system. When T = Z and T = Kg we get the lattice and g-discrete counterparts 
of the AKNS hierarchy where the first nontrivial flow is (1371) . Besides, for T = Z the system 
f l37|) . together with its bi-Hamiltonian structure, is equivalent to the system considered in [9]. 



A-differential Kaup-Broer, k = 1. Consider the following Lax operator with the con- 
straint dMD (a = 1) 

L = (1 + /iw — jj^w)5 + + 5~^w. 
The flrst and the second flows are 

= (1 + yui; — ^^w)Av — fiA\w + fivw — ^w'^), ^^^^ 
Wt^ = —At (u7 + fivw — jJ^w"^^ 

and 

= uA {iP' + 2uw + uAv + fiA\uw)) — fiA^ (2uvw + fxuwAv + uA\uw)) ^^^^ 
= —At {2uvw + jjiuwAv + 'uAt('uw)) , 



where 



The differentials are given in the form 



dH = 6 ^— + ^ + yu— . 

ov ow ov 



Thus, the hnear Poisson tensor flMj) is 

-A^u 
The recursion operator has the form [IS] 

w + uA + R fiv — fiw + (2 + fiA'^)u — Rfj, 



where 
Hence 



w — A'^uwA ^ A'^u + V — fiwA'^uwA ^fiu ^ 
R = uAvA'^rr^ - fiAUlwA^^u'^. 



_ ^ _ f TTyv uAv + uAuA + fiuwA — ^A'^uw\ 

'^^ \^— t;A''''S — A^uA^u + uwAfi — A^fiuw uwA — A'^uw J 

where 

TTyy = uAfiv — /if A^M + uAu — uA^u + uAuAfi — fiA^uA^u + fiuwA^ — fiA^ ^uw. 
The Hamiltonians are 

Hq = I w Ax 
Jj 

Hi = J ^cw — 2''^^^^ 

H2 = J + v'^w + {w + fivw — fi'^w'^)Av — -fi^w^^ Ax 



When T = M the above construction recovers the field Kaup-Broer hierarchy with its bi- 
Hamiltonian structure [H] . As previously, the Kaup-Broer system is given only by the second 
flow ( HTl) . In the lattice case, i.e. of T = Z, the above bi-Hamiltonian hierarchy is equivalent 
to the relativistic Toda hierarchy considered in [9] and, (HOl) is equivalent to the relativistic 
Toda system. 



8 Conclusions 

We have presented a unified theory of the construction of the bi-Hamiltonian nonlinear evo- 
lution hierarchies such as field, lattice and g-discrete soliton hierarchies. Actually, we took 
advantage from the theory of time scales. Therefore, one can also consider the construction of 
soliton systems with spatial variable belonging to the spaces being partially continuous and 
discrete. This might be interesting from the point of view of applications. There are also 



other approaches generahzing and unifying theory of sohton systems, presented in [22] and 

On the other hand, making use of the regular-discrete time scales only, the theory from the 
article can be considered as a discretization scheme of field soliton systems. In some special 
cases, introducing appropriately deformation parameter to some regular-discrete time scales, 
one can consider the quasi-classical limit of discrete soliton systems yielding dispersive field 
soliton equations. In particular for T = KL the quasi-classical limit is given by — and for 
T = Kg by g ^ 1, see [I5]. 



Acknowledgement 

This work was partially supported by the Scientific and Technical Research Council of Turkey 
and MNiSW research grant no. N N202 404933. B.Sz. was supported by the European 
Community under a Marie Curie Intra-European Fellowship, contract no. PIEF-GA-2008- 
221624. 



Appendix 

We verify the Lemma 14.11 by considering the positive and negative cases of n separately by 
the use of induction. Note that E — yuA = — fiA^ = 1. 

Let n ^ 0. Assume that (fT7|) holds for positive n. Then 
{E - /iA)"+i = {E- fiAy\E - fiA) = {E- fiAy'E - fi{E - fiA^A 

n n n+1 n+1 

= + Y.^-f^)'^'S^A = Y.{-fi)'S^E + Y,i-f^)'S^_,A 

k=0 k=0 k=0 k=0 

n+1 n+1 

k=0 k=0 

where we used the fact that S^^j^i = S^i = and the recurrence relation ( ITSl) . 

Let n < 0. First we show ( |T5i) for n = —1. Thus, using the recursive substitution we 
deduce 

{E - /iA)"i = {E-^ - i^iA^) {E - /iA)-i = E^^ - fi{E - fiA)-^A^ 

= E~^-I2 {E~^ - i2{E - fiAy^A^) At = E~^ - fiE^^A^ + ^^(^ _ i^A^A^'' 
= E-' - fiE-'A^ + fi^E-'A^' - /i^E-i At' + . . . 

oo cx> 
k=0 k=0 



Assume that ffTTl) holds for negative n. Then 

oo oo 

{E - /iA)"-i = {E- fiATiE - /iA)-i = Y.{-fi)'Sj: Y.{-f,yE''A^' 

fc=0 i=0 

oo oo oo fc 

fc=0 i=0 fc=0 i=0 

oo fc oo 



fc=0 j=0 fc=0 



where we used fll7l) for n = — 1 and the recurrence relation (fTBj) . Hence fll7l) holds for n — 1, 
which finishes the proof. 
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